Enhancing or suppressing spin Hall effect of light in layered nanostructures 
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The spin Hall effect (SHE) of light in layered nanostructures is investigated theoretically in this 
paper. A general propagation model describing the spin-dependent transverse splitting of wave 
packet in the SHE of light is established from the viewpoint of classical electrodynamics. We show 
that the transverse displacement of wave-packet centroid can be tuned to either a negative or a 
positive value, or even zero, by just adjusting the structure parameters, suggesting that the SHE of 
light in layered nanostructures can be enhanced or suppressed in a desired way. The inherent physics 
behind this interesting phenomenon is found to be attributed to the optical Fabry-Perot resonance. 
We believe that these findings will open the possibility for developing new nano-photonic devices. 

PACS numbers: 42.25.-p, 42.79.-e, 41.20.Jb 



I. INTRODUCTION 

Spin Hall effect (SHE) is a transport phenomenon, in 
which an applied field on the spin particles leads to a 
spin-dependent displacement perpendicular to the elec- 
tric field direction The SHE of light can be re- 
garded as a direct optical analogy of SHE in electronic 
system where the spin electrons and electric potential are 
replaced by spin photons and refractive index gradient, 
respectively [J-Q . The SHE of light is sometimes referred 
to as the Fedorov-Imbert effect, which was predicted the- 
oretically by Fedorov Q, and experimentally confirmed 
by Imbert [8| . The spin-dependent transverse shift in the 
SHE of light is generally believed as a result of an effec- 
tive spin-orbital interaction, which describes the mutual 
influence of the spin (polarization) and trajectory of the 
light beam [§]. 

Recently, the SHE of light has been extensively inves- 
tigated in different physical systems. In a static gravi- 
tational field, the photon Hamiltonian shows a new kind 
of helicity-torsion coupling, resulting in a novel birefrin- 
gence phenomenon: photons with distinct helicity follow 
different geodesies [10|]. In optical systems, the SHE of 
light is observed directly in the glass cylinder and its 
fundamental origin is related to the dynamical action of 
the topological Berry-phase monopole in the evolution 
of light [ll|. The SHE of light can also be observed in 
scattering from dielectric spheres [Hj]. In particular, a 
giant SHE of light can be produced by subwavelength 
displacements of a nanoparticle fl3j . Even in free space, 
the SHE of light can be observed on the direction tilted 
with respect to beam propagation axis [Tij . In plasmonic 
systems, a spin-dependent splitting of the focal spot of a 
plasmonic focusing lens was demonstrated and explained 
in terms of a geometric phase flEl ]. In semiconductor 
physics, the SHE of light has been observed in silicon via 
free-carrier absorption. The interesting result suggests 
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that the SHE of light has the potential of probing spatial 
distributions of electron spin states [n| . 



The SHE may offer an effective way to manipulate the 
spin particles, and open a promising way to some po- 
tential applications, such as in dense data storage, ultra- 
fast information processing, and even quantum comput- 
ing [l7l - fl{| . The generation, manipulation, and detection 
of spin-polarized electrons in semiconductors and nanos- 
tructures define the main challenges of spin-based elec- 
tronics. Similar challenges also exist in spin-based pho- 
tonics. The SHE of light may open new opportunities for 
manipulating photon spin and developing new generation 
of all-optical devices as counterpart of recently presented 
spintronics devices. In this paper, we will study the SHE 
of light in layered nanostructures in which the refractive 
indices of their constituent materials vary between high- 
index regions and low-index regions. Such an environ- 
ment presents to photons as an analogy of semiconduc- 
tor presenting potential to electrons |20|], thus present- 
ing some imaginable interesting properties of the SHE of 
light. 



The paper is organized as follows. First, we want to 
establish a three-dimensional propagation model to de- 
scribe the SHE of light in layered nanostructure. The 
Fresnel coefficients are no longer real in the layered 
nanostructures, so it is necessary for us to obtain a 
more general expression. Next, we attempt to reveal 
what roles the Fresnel reflection and transmission co- 
efficients play in the SHE of light. We find that the 
Fresnel coefficients present sine-like oscillations and the 
spin-dependent splitting of wave-packet centroid signifi- 
cantly depends on their ratio. Finally, we want to explore 
the secret underlying this interesting phenomenon. The 
result shows that the SHE of light can be readily modu- 
lated, i.e., enhanced or suppressed, via tuning the optical 
resonance in layered nanostructures. 
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II. THREE-DIMENSIONAL BEAM 
PROPAGATION MODEL 

Figure [T] illustrates the beam reflection and refraction 
in the layered nanostructure. The z axis of the lab- 
oratory Cartesian frame (x,y,z) is normal to the in- 
terfaces of the layered structure. We use the coordi- 
nate frames (x a ,y a , z a ) for central wave vector, where 
a = i,r,t denotes incident, reflection, and transmission, 
respectively. We apply the angular spectrum method to 
derive an expression for a three-dimensional beam prop- 
agation model. Hence, we use local Cartesian frames 
(X a ,Y a , Z a ) to describe an arbitrary angular spectrum. 
The electric field of the ath beam can be solved by em- 
ploying the Fourier transformations. The complex am- 
plitude for the ath beam can be conveniently expressed 
as [13 

E a (x a , y a , z a ) — J dk ax dk a y^ a {k aXl & a y) 
x exp[i(k 



where k az = Jk%- {kl x + k% y ) and E a (fc QX , k ay ) is the 

angular spectrum. The approximate paraxial expression 
for the field in Eq. (fT]) can be obtained by the expansion 
of the square root of k az to the first order [22j, which 
yields 



E a — Gxp(ik a z a ) I dk ax dk ay Ei a {k ax ,k a y) 



x cxp 



2k a 



■(2) 



In general, an arbitrary linear polarization can be de- 
composed into horizontal and vertical components. In 
the spin basis set, the angular spectrum can be written 
as: 
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Ef = — (E 
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-^i(E ( - - E, 



(3) 



(4) 



Here, H and V represent horizontal and vertical polar- 
izations, respectively. The positive and negative signs 
denote the left and ri ght circularly polarized (spin) com- 
ponents, respectively [23|]. The monochromatic Gaussian 
beam can be formulated as a localized wave packet whose 
spectrum is arbitrarily narrow, and can be written as 



Ei± = —7={eix ± ie iy)—7== exp 
v 2 v 27r 



™U4 



kfy) 



, (5) 



where wq is the beam waist. After the angular spectrum 
is known, we can obtain the field characteristics for the 
ath beam. 

To accurately describe the SHE of light in layered 
nanostructure, it is need to determine the reflection and 




FIG. 1: (color online) Schematic illustrating the reflection and 
refraction of central and local wave vectors at (m+ l)-layered 
nanostructure. x a y a z a (a = i,r,t) are reference frames for 
central wave vector of incident, reflection, and transmission, 
respectively, oxyz and OXYZ are the interface reference 
frames for central and local wave vectors, respectively. 



transmission of arbitrary wave- vector components, which 
can be solved by 2 x 2 transmission matrix [24j ]: 

M = ToiPiTi2P2---Pm-2T m -2,m-lPm-lT m ~l, m , (6) 

where 



t, 



1 T m -\. 
Tm — l.m 1 



(7) 



is the transformation matrix from (to — l)-th to m-th 
layer, and 



Pn 



exp[ik mz d m 




exp(- 







(8) 



is the transmission matrix for m-th layer. Here, d m is 
the thickness of TO-th layer, r m _i iin and t m _i, m are re- 
flection and transmission coefficients from (m — l)-th to 
TO-th layer, respectively. For an arbitrary wave-vector 
component, the Fresnel coefficients of the layered nanos- 
tructures can be written as 



1 



H21 



P>8 



(9) 



where p and s denote parallel and perpendicular polar- 
izations, respectively. By making use of Taylor series 
expansion, the Fresnel coefficients can be expanded as a 
polynomial of ki X . We obtain a sufficiently good approx- 
imation when the Taylor series are confined to the zero 
order. 

From the boundary condition, we obtain k rx = — ki X 



and 



k — 

t^ry — 



After a series of calculations of the 



reflected angular spectrum given in the Appendix [XJ 
Eq. © together with Eqs. © and (|a7|) provides the 
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paraxial expression of the reflected field: 



E 



H 



ZR 



x exp 



7rt«o 2fl + iz r 
k x* + (y r ±5?) 2 ~ [ 



exp(ik r z r ) 



zr + iz r 



(10) 



E 



ZR 



x exp 



exp(ifc r z r ) 
,2 , /„, i jtV<\2- 



fc ^ + (y r ±<5;)^ 



(11) 



where = k wl/2 is the Rayleigh lengths, = (f + 
r s /r p ) cot 9i/k and <^ = (f + r p /r s ) cot ^/fc . 

We next consider the transmitted field. From the 
Snell's law under the paraxial approximation, we obtain 
ktx = ki X /rj and kt y = h y . Substituting Eqs. © and 
(|A10p into Eq. ©, we obtain the transmitted field: 



E 



H 



t p (e tx ± ie ty ) 



z Ry exp(ik t z t ) 



x exp 



™0 \/(z Rx + iz t )(z Ry + iz t ) 
n m k ( xj (y t T 5?) 2 

zr x + iz t ZR V + izt 



1-2) 



E 



^it s (e tx ± ie ty ) z Ry exp(ik t z t ) 



t± 



^ w o \J{.zrx + izt)(z Ry + iz t ) 



x exp 



n m k() 



zr x + iz t ZRy + iz t 



13) 



Here. 



S? 



(?] — t s /t p )cot9i/k and 6} 



(v 



tp/ts) cot8i/ko. The interesting point we want to stress 
is that there are two different Rayleigh lengths, zr x = 
nrnV 2 koWQ / 2 and zr v = n m k§w\l2, characterizing the 
spreading of the beam in the direction of x and y axes, 
respectively (25|. Note that the Fresnel coefficients are 
no longer real in the layered nanostructure. Hence, we 
should extend the previous expression of transverse dis- 
placement [6j to a more general situation. 



III. SPIN HALL EFFECT OF LIGHT 

It is well known that the SHE of light manifests itself 
as polarization-dependent transverse splitting. To reveal 
the SHE of light, we now determine the transverse dis- 
placements of field centroid. The time-averaged linear 
momentum density associated with the electromagnetic 
field can be shown to be 126] 



p„(r) = ^Re[E a (r) x H^r) 



(14) 



where the magnetic field can be obtained by H a = 
— ik~^ x E a . The intensity distribution of wave packet 
is closely linked to the longitudinal momentum currents 

lyKai Va,-* %a) ^ Pa ' ^ az . 



At any given plane z a = const., uui, .mnor 
placement of wave-packet centroid compared 
geometrical-optics prediction is given by 



the transverse dis- 
to the 



J / y a I(x a ,y ai z a )dx a dy a 
J J I(x a ,y a ,z a )dx a dy a 



(15) 



Note that the transverse displacement can be divided into 
z a -dependent and z a - independent terms. We here con- 
centrate our attention on the z Q -indepcndent transverse 
displacements. 

We first consider the spin-dependent transverse dis- 
placement of the reflected field. After substituting the 
reflected field Eqs. (dUJ) and (HU) into Eq. ([T5j) , we obtain 
the transverse spatial displacements as 



4" 



i(ip s - ip p )] COt( 



Ttt-II + \r p \/\r s \cos(ip p - (p s )] cotf 
Ztt 



(16) 



(17) 



where r P:S = |r PjS | exp(«(^ PiS ). Note that these expres- 
sions are slightly different from the previous work [6J, [27J , 
since the Fresnel reflection coefficients are no longer real 
in our model. 

We next consider the spin-dependent transverse dis- 
placements of the transmitted field. After substituting 
the transmitted field Eqs. (H2]) and (13]) into Eq. ((T5|), we 
have 



Ayf ± 



AyY ± 



A . 

Z7T 



\tp\ cos( 



±7^^- \tp\/\t s \cOs((t>p 



^p)]cot0i, (18) 



^)]cot0i, (19) 



where i PjS = \t pa \ exp(j0 PjS ). For an arbitrary linearly 
polarized incident beam, the calculation of the transverse 
displacements for the reflected and transmitted field is 
given in the Appendix [Bj 

For left and right circularly polarized components, the 
eigenvalues of the transverse displacement are the same 
in magnitude but opposite in directions. Under the limit 
of m = 1 (air-glass interface), the above expression coin- 
cides well with the early results 0. Our scheme shows 
that the SHE of light can be explained from the viewpoint 
of classic electrodynamics. For incidence angles greater 
than the critical angle of total internal reflection, most of 
photons are reflected, and part of them tunnel through 
the layered structure [28[ . Hereafter, we only concentrate 
our attention on the SHE of light in the transmission case. 

As an example, a five-layered nanostructure composed 
of two prisms (Pi and P2), two films (MgF2), and an air 
gap [Fig. |2Ja)] is chosen to illustrate the SHE of light. 
We consider two kind of systems: (i) A symmetric sys- 
tem with two BK7 prisms (uq — — 1.515 at 633 nm) 
filmed by MgF2 [ri\ = 713 = 1.377 at 633 nm) and the 
middle layer is air gap (r%2 — 1); (ii) Replacing P2 by 
an S-LAH79 prism (n.4 = 1.996 at 633 nm) to form an 
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FIG. 2: (Color online) The Fresnel transmission coefficients 
in a five-layered nanostructure. (a) Schematic illustrating 
the layered nanostructure: The two prisms (Pi and P2) both 
filmed by a low refractive index layer with n\ = 713 = 1.377 
(MgF2 at 633 nm), di = d3 = 80nm, and the thickness of the 
air gap (tfe) can be modulated. The incident angle is chosen 
as 9i — 7r/6. The Fresnel transmission coefficients t p and t s 
versus the thickness of the air gap d,2 for the the symmetric 
and asymmetric system: (b) no = ru = 1.515 (BK7 at 633 
nm) and (d) n = 1.515 (BK7 at 633 nm) and n 4 = 1.996 (S- 
LAH79 at 633 nm). (c) and (e) are the phase of the Fresnel 
transmission coefficients for the two systems, respectively. 



asymmetric system. The MgF2 layers with 80 nm thick- 
ness were prepared on Pi and P2. For a given incident 
angle Qi = tt/6, t s in the two systems both behave sine- 
like oscillations versus to the thickness of air gap (cfe) due 
to the optical Febry-Perot resonance with multi-resonant 
peaks for different cfo, while t p is nearly unchanged and 
insensitive to di [Fig. Hfb) and^d)]. It is obvious that 



|i P I cos( 



b p ) and |i p |/|i s | cos( 



in Eqs. (TT 



and (1T91) determine the magnitude of the transverse dis- 
placements (Ayt±) of wave-packet centroid since other 
quantities in the two equations are constant. Actually, 
the cj> p is nearly equal to <fi s versus (L2 for both the sym- 
metric [Fig. EI[c)] and asymmetric system [Fig. Hfe)], 
which means cos(0 p — 4>s) ~ 1 and the magnitude of 
Ay t ± only depends on |t s |/|£ p | or |t p |/%|. 

The spin-dependent splitting in SHE of light is 
schematically shown in Fig. (2 a )- From Eqs. (TT5)) 
and (|19l) . we know that the transverse displacements of 
H and V polarizations would have just opposite tendency 
versus di which can also be seen from Fig.[3Jb)-[3][e). For 
a fixed incident angle of Qi = tt/6, the transverse dis- 
placements present sine- like oscillations since the value of 
|ts|/|t P | or |t p |/|t s | is periodic due to the Fabry- Perot res- 
onance in the layered structures. In the symmetric case, 
the transverse displacements present a sine-like oscilla- 
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FIG. 3: (Color online) (a) Schematic illustrating the spin- 
dependent transverse splitting. The transverse displacements 
versus the thickness of the air gap c?2 in the situation of trans- 
mission with incident angle Qi — tt/6. (b) and (c) show the 
transverse displacements of H and V components, respec- 
tively, for the symmetric system, (d) and (e) also represent 
the transverse displacements but for the asymmetric system. 



tion in the range of zero-positive-zero or zero-negative- 
zero for a certain polarization component [Fig.[3jb)-[3jc)]. 
In the asymmetric case, for ii-polarization, the trans- 
verse displacements can exhibit a positive and negative 
transverse displacement for left- and right-polarized com- 
ponents, respectively [Fig. GUd)]. It is interesting to note 
that the transverse displacements exhibits a sine-like os- 
cillation in the range of negative-zero-positive values for 
V^-polarization, which can be modulated via tuning the 
thickness of the air gap [Fig. Efe)]. It indicates that the 
SHE of light can be greatly enhanced or suppressed, or 
even completely eliminated. 

From the above analysis, we know that the transverse 
displacements are related to the ratio between the Fres- 
nel transmission coefficients |t p | and \t s \, whose depen- 
dence on the thickness of the air gap are periodic due 
to the Fabry-Perot resonance in the layered nanostruc- 
tures. Hence, we can expect that a nanostructure with 
a large ratio of |i a |/|tp| or |i p |/|i s | would extremely en- 
hance the SHE of light. On the contrary, a small ratio 
of |t s |/|ip| or | t p |/| t s | would greatly suppress the SHE of 
light. In fact, the phase difference <fr p — (j> s of the Fresnel 
transmission coefficients may also be used to modulate 
the SHE of light since it can change the sign of Ayt±. 
The nanostructures with a metamaterial layer whose re- 
fractive index can be tailored arbitrarily can be a good 
candidate to support this prediction [25| . Under the con- 
dition of |i s |/|tj,| cos(</> s — 4> p ) — tj (H polarization) and 
l*p|/l*«l cos((p p — 4> s ) — rj (V polarization), the SHE of 
light can be suppressed completely. 

It should be mentioned that the spatial separation 
of the spin components can also be tuned continuously 
by varying the incident angle in a single air-glass inter- 
face Q. However, the refracted angle changes and the 
transmission coefficients decrease accordingly as the in- 
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cident angle increases. Hence, it is disadvantage for po- 
tential application to nano-photonic devices. As shown 
in above, the wave packet in our scheme is incident at 
a fixed angle, and the transverse displacements can be 
tuned to either a negative or a positive value, or even 
zero, by just adjusting the structure parameters. Mean- 
while, the wave packet in the layered nanostructure ex- 
hibit much higher transmission coefficients than in the 
single air-glass interface. Hence, the layered nanostruc- 
tures provide more flexibility for modulating the SHE of 
light. We will search for suitable nanostructures to ma- 
nipulate it in the future. It is expected that the SHE of 
light in layered nanostructures will be useful for design- 
ing very fast optical switches, for example, by replacing 
the air gap by material whose refractive index can be 
tuned by a electric field. 



Then, we transform the electric field from the reference 
frame xyz around the y axis by an angle ki y /(ko sin Oi) to 
the frame XYZ, and the correspondingly matrix is given 

by 



m xyz -+xYZ 



1 

kiy 



ko sin Oi 

o 



fep sin ( 

l 
o 



(A2) 



where fco is the wave number in vacuum. Finally, we 
transform the electric field from the reference frame 
XYZ around the y axis by an angle —Oi to the frame 



XiYiZi, and the matrix can be written as 



niXYZ-t-XiYiZi — 



cos Oi sin 6 
1 
— sin Oi cos ( 



(A3) 



IV. CONCLUSIONS 

In conclusion, we have revealed a tunable SHE of light 
in layered nanostructures. From the viewpoint of classi- 
cal electrodynamics, we have established a general prop- 
agation model to describe the spin-dependent transverse 
splitting of wave packet in the SHE of light. By modu- 
lating the structure parameters, the transverse displace- 
ments exhibit tunable values ranging from negative to 
positive, including zero, which means that the SHE of 
light can be greatly enhanced or suppressed, or com- 
pletely eliminated. We have shown that the physical 
mechanism underlying this intriguing phenomenon is the 
optical Fabry-Perot resonance in the layered nanostruc- 
ture. These findings provide a pathway for modulating 
the SHE of light, and thereby open the possibility for 
developing new nano-photonic devices. 
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Appendix A: Calculation of the reflected and 
transmitted angular spectra 

In this appendix we give a detailed calculation of the 
reflected and transmitted angular spectra. From the cen- 
tral frame XiyiZi to the local frame XiYiZi, the following 
three steps should be carried out. First, we transform 
the electric field from the reference frame Xiy t Zi around 
the y axis by the incident angle 0i to the frame xyz: 



E 



xyz 



>xyzE XiyiZi , where 



^XiyiZi—txyz 



cos 0i — sin 

1 
sin 0i cos 0i 



(Al) 



Thus, the 
tral frame 
can be 



rotation 
XiyiZi to 
written 



matrix from the cen- 
the local frame XiYiZi 



as 



M, 



XiyiZi 



> XiYiZi 



THxYZ^x i Y i z i 'm X y Z ^ t xYzm Xi y iZi ^ tX y Z , and we have 



Al 



XiyiZi^t XiYiZi 



1 

kiy cot 0i 
ko 



k iy cot 6 
ko 
1 

kiy 

kg 





k 
1 



(A4) 



For an arbitrary wave vector, the reflected field is de- 



E XiYiZi, where r p and 



termined by E Xr Y r z r = r PyS 
are the Fresnel reflection coefficients. The reflected field 
should be transformed from X r Y r Z r to x r y r z r . Fol- 
lowing the similar procedure, the reflected field can be 
obtained by carrying out three steps of transformation: 



E. 



f U f % f 



Al 



E Xr Y r z r where 



M 



X r Y r Z r 



k ry COt 0i 

ko 



k ry cot 0i 
ko 
1 



(A5) 



Here, only the two-dimensional rotation matrices is taken 
into account, since the longitudinal component of elec- 
tric field can be obtained from the divergence equation 
E az kaz = -(Eaxkax + E ay k ay ). The reflection matrix 
can be written as 



Mr = M Xr Y r Z r - 











AL 



xtytZi^XiYiZi 



(A6) 



The reflected angular spectrum is related to the bound- 
ary distribution of the electric field by means of the re- 
lation E r (k rxi k ry ) = MnEi(ki Xl kiy), and we have 



k ry (r p -\-r s ) cot $i 



kry (r p -\-r s ) cot 0i 
ko 



Ef 

eY 



A7) 



We proceed to consider the transmitted field. Follow- 
ing the similar procedure, we obtain the transform matrix 
from X t Y t Z t to x t y t z t as 



M 



XtY t Zt->x t y t z t — 



k ty cos t 
ko sin 0i 



k ty cos 6> t 
ko sin 0i 
1 



(A8) 
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where 9 t is the transmitted angle. For an arbitrary wave 
vector, the transmitted field is determined by Ex t Y t z t = 
tp^ExiYiZi, where t p and t s are the Fresnel transmis- 
sion coefficients. Hence, the transmission matrix can be 
written as 



M T - M Xt Y t Z t ->x t y t 2 



t p 
t. 



M 



XiyiZi^XiYiZi 



(A9) 



The transmitted angular spectrum is related to the 
boundary distribution of the electric field by means of 
the relation E t (k tx ,k ty ) = MTEi{ki X , ki V ), and can be 
written as 







. E Y _ 





tp 

kty (rjtp—ts) cot 6i 
feo 



kty(t p —r)t s ) COt Oj 
t. 



eY [ 



where rj — cos6 t / cos Qi. 



Appendix B: Transverse displacements for arbitrary 
linear polarization 

For an arbitrary linearly polarized beam, the trans- 
verse displacements of the reflected field are given by 

Ay r± = cos 2 ■y r AVr± + sin2 7rAy^ ± , (Bl) 

where j r is the reflected polarization angle. In the frame 
of classical electrodynamics, the reflection polarization 



angle is determined by: 



cos7 r 



cos7iRe[r p ] 



sin J r 



cos 2 7iRe[r p ] 2 + sin 2 ^Ke [r s ] 2 

sin7 t Rc[r s ] 

^ cos 2 7iRe[r p ] 2 + sin 2 7iRe[r s ] 2 



(B2) 



(B3) 



Here, 7, is the incident polarization angle. For an arbi- 
trary linearly polarized wave-packet, the transverse dis- 
placements of the transmitted field are given by 

Ay t ± = cos 2 7*Ai/& + sin 2 7t AyY ± , (B4) 

where the transmission polarization angle 7* determined 

by 



cos 7* = 



C0S7iRe[t p 



sin7 t 



\J cos 2 7iRe[tj,] 2 + sin 2 jiRe[t s } 2 

sin7iRe[t s ] 
cos 2 jiRe[tp} 2 + sin 2 7iRe[i s ] 2 



(B5) 



(B6) 
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